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Abstract. Modifying the standard approaches to nonperturbative QCD based on Borel-transformed 
dispersive sum rules by allowing the effective continuum thresholds required for the implementation 
of quark-hadron duality to depend on the Borel parameters and on any relevant momentum promises 
to entail higher accuracy and reliable error estimates for the extracted predictions of hadron features. 
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INCENTIVE: IMPROVEMENT OF QCD SUM RULES [1-14] 

Within the method of QCD sum rules, a la Shifman-Vainshtein-Zakharov, the concept of 
^ — ■ quark-hadron duality is usually implemented by assuming that above specific continuum 

\Q '. thresholds the contributions to suitably defined correlators of interpolating currents at the 

level of the QCD degrees of freedom equal those at the level of the hadronic bound states. 
Our approach [1-14] seeks to quantify the uncertainty induced by such an approximation 
and to improve the accuracy of all predictions by allowing our threshold to depend on the 
involved momenta and parameters introduced upon application of Borel transformations. 
The Borel transformation to a new variable, the Borel parameter (called T or x hereafter), 
serves to remove subtraction terms and suppress excitation and continuum contributions. 

Any such idea is best tested first in a situation in which the outcome for all bound- state 
characteristics one is interested in is known exactly. So let's study a quantum-mechanical 
model defined by a nonrelativistic Hamiltonian H with a harmonic-oscillator interaction: 
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This model is exactly solvable: ground-state energy, decay constant, and form factor read 
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We inspect the analogues of a few Af-point correlation functions of interpolating currents. 
The Borelized polarization function or 2-point vacuum-vacuum correlator [1-5, 8] reads 



n(T) = (x f = 0| exp(-// T) | Xi = 0) ™ 
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The double-Borelized 3-point correlator of some current operators 7(q) [6, 9] is given by 

r(72,r 1 , g ) = (x f =o| e xp(- J ffr 2 )/(q)exp(-/fr 1 )|xi = o) 

R ge xp[-E g (T i+ T 2 )]F g (q); 
for equal Borel parameters (or Euclidean times) T\=Ti = \ this expression simplifies to 



r(r,<7) =°n(r)exp 
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R g exp(-£ g T) F g {q) = T g (T, q) 



The Borelized vacuum-hadron amplitude of the T-product of 2 'quark' currents is [7, 12] 



A(r,4) = (x = 0|exp(-tf7V(q)|v/ g > 



%exp(-E g T)F g (q) = A g (T,q) 



The admissible working ranges of Borel parameter values, or 'Borel windows' , are found 
by requiring that the ground state contributes sizeably to the average energies (see Fig. 1) 

E T( T, q) = -^iogr ( r, 9 ) = | m co,h( m r) + 4m[l+ J o 2 sh(mr)1 , 
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E A (T,q) = - \ og A(T,q) = - CO + |^exp(-2 COT) , 
ol 2 2m 

and that in the sum rules higher-order nonperturbative effects are sufficiently suppressed. 
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FIGURE 1. Relative ground-state contributions to Er(T,q) andr(r,<7) (first row), respectively Ea (T, q) 
and A(T,q) (second row), for several dimensionless q = qj y/mo5; vertical lines delimit the Borel windows. 



DUALITY EFFECTIVE CONTINUUM THRESHOLDS 



The derivation of any QCD sum rule proceeds along several more or less canonical steps: 

• Perform the operator product expansion (OPE) of the T-product of involved fields in 
terms of local operators yielding so-called vacuum condensates. In quantum theory, 
an OPE corresponds to a series expansion in powers of the Borel parameter [1-3, 8]. 

• Represent the (perturbative) correlator as dispersion integral over a spectral density. 

• Invoke quark-hadron duality, by cancelling the contributions to the correlator above 
its effective continuum threshold against its hadronic 'continuum' : we argue that all 
thresholds must depend on the Borel parameter and, where applicable, on momenta. 

Even after application of the quark-hadron duality approximation the resulting sum rules 
can be rendered rigorous by introducing the notion of dual correlators. These QCD-level 
correlators are defined to be exactly dual to the relevant lowest, or ground- state, hadronic 
contributions. This becomes feasible if allowing the exact effective continuum thresholds 
ZeS (T, q) to depend on the Borel parameter (and, if relevant, on the involved momentum): 

n d uai(7\z? ff (t)) =Jdzexp(-zT) Po (z) +n power (r) s = n g (r) =i? g ex P (- J E g r) , 
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r dnai (T,q,zl s (T,q)) =Jdzi J dzi exp ( - Zl * Zl T j A (zi ,Z2,q) + T powei (T,q) 
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= r g (r,4) = R g exp(-E g T)F g (q) , 

Adwd(T,q,zts(T,q)) = J dzexpf-zT) p A (z,T,q) =A g (T,q) = v /7? g exp(- J E g r) F g (q) . 



Our exact knowledge of the quantum-mechanical ground-state mass, decay constant and 
form factor renders possible to determine the exact effective continuum thresholds for all 
dual correlators: Fig. 2 shows that these thresholds are definitely far from being constant. 
However, in order to extract from rd ua i(7\ q, z^ s (T, q) ) or Ad ua i(7\ q, Z^ & (T, q) ) the elastic 
form factor F g (q) by exploiting only the knowledge of E g and R g we simultaneously have 
to pin down the actual behaviour of the effective continuum thresholds. We thus adopt for 

our thresholds the simple polynomial ansatz [10-12] z^(T,q) = T!j=ozf\q) (ooT)-> and 
obtain any set of coefficients z^ (q) by evaluating the dual energies, defined according to 

^duai^?) = -^iogr dual (r,4,z|rf f (r,4)) , 

ELl(T,q) = --^\ogA dnal (T,q,4 s (T,q)) , 

at several Borel-parameter points 7}, 1 < i < N, within the Borel window and minimizing 
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then the squared difference X 2 = jj I4L 1 £ dua i (Ti,q)-E g between £ dual ( 7] , q) and E g . 
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FIGURE 2. Exact effective continuum thresholds Zg Xact (T,q) (first row) and z^ xact (T,q) (second row) for 
r dml (T : q,z r eS (T,q)) <mdA dual (T,q,4 n (T,q)): r^^O) =n(r!+r 2 ) yields zl xact (T,q = 0) = Z? xact (T). 

Figure 3 depicts the Borel behaviour of our resulting dual energies and dual form factors. 
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FIGURE 3. Fitted dual energy £duai (T, q) (left) and dual elastic form factor Fd ua i (T, q) (right) for T(T, q) 
(upper row) and A(T,q) (lower row) at dimensionless q = q/y/mco — 2. These plots reveal that the intrinsic 
uncertainties of the sum-rule method cannot be reliably estimated if assuming a constant threshold (n = 0). 



BRIDGING THE GAP: QCD [13] < — > POTENTIAL MODELS [14] 



Beyond doubt, our ultimate target must be to apply our insights on the sum-rule approach 
gained in the course of the above quantum-theoretical exercise to QCD in order to extract 
characteristics of actual hadrons such as decay constants and form factors. Let us start by 
studying the Borelized correlator of 2 pseudoscalar quark currents 75 = (m^ + m u ) u i 75 b: 
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n(T) =ydsexp(-ST)p p e r t(s,M)+n p ower(T, i u) 
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The perturbative spectral density p pe rt(s, M) ma y be derived from QCD in form of a series 
expansion in powers of the strong coupling oc s , at an appropriate renormalization scale /1: 

f \ / \ aJu) , . ocl(u) , \ 

Ppert(j,M) = Po(s,M) + ^T^Pl( s >i U ) + ^pP2(^M) + -- - ■ 
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FIGURE 4. Dual masses and decay constants extracted from QCD (left) and our potential model (right). 
Duality permits dual correlators exactly counterbalancing all ground-state contributions: 
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The decay constant /b of the charged B meson is defined by (mb + m u ) (0| 7*5 |B) = fsM^. 
To parallel the quantum-theoretical case we introduce dual B mass and decay constant by 

M dual( T ) = logn dual (T, Sefif(T)) , /dual( T ) = ^ B 4ex P( M B n dual (T, s eff (T)) . 

We make a polynomial ansatz for the effective continuum threshold s[^ ( t) = L"=o s ^ tJ 

and get the coefficients by minimizing X 2 = jj E;li [^duai( T — ^b] 2 f° r T *<N m me 
Borel window [13]. Figure 4 confronts the QCD answer for the B meson with that from a 
quantum-theoretical model with harmonic-oscillator plus Coulomb interaction potential. 

OBSERVATIONS, RESULTS, FINDINGS, AND CONCLUSIONS 

1 . Our exact effective continuum thresholds do depend on the Borel parameter and any 
relevant momentum, and they are not universal but vary with the studied correlators. 

2. The algorithm proposed here both improves significantly the accuracy of traditional 
sum-rule predictions and provides reliable estimates of their intrinsic uncertainties. 

3. The striking (and surprisingly even quantitative) similarity of our hadron-parameter 
extraction procedures in potential models and in QCD calls for application in QCD. 
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